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Abstract: We investigate the dependence of the conductivity and the entanglement en¬ 
tropy on the space-time dimensionality d in two holographic superconductors: one dual 
to a quantum critical point with spontaneous symmetry breaking, and the other modeled 
by a charged scalar that condenses at a sufficiently low temperature in the presence of 
a Maxwell field. In both cases the gravity background is asymptotically Anti de Sitter 
(AdS). In the large d limit we obtain explicit analytical results for the conductivity at zero 
temperature and the entanglement entropy by a 1 /d expansion. We show that the entan¬ 
glement entropy is always smaller in the broken phase. As dimensionality increases, the 
entanglement entropy decreases, the coherence peak in the conductivity becomes narrower 
and the ratio between the energy gap and the critical temperature decreases. These results 
suggest that the condensate interactions become weaker in high spatial dimensions. 
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1 Introduction 

It is a well known fact in condensed matter and statistical physics that the dynamics of 
many systems simplifies drastically in the limit of large space dimensions d — 1 [1-9]. 
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Analytical results are typically obtained for d —> oo [1, 2] and, in some cases, it is also 
possible to compute explicitly small corrections [7] due to a large but finite dimensionality 
by a 1/d expansion. A typical example is the Hubbard model in the strong coupling 
region where in the large d limit the problem maps onto a mean-field quantum impurity 
model that is solved self consistently. Meaningful results are only obtained [1] after the 
kinetic energy is properly rescaled so that the trivial non-interacting limit is avoided for 
d —> oo. The application of these ideas to the Hubbard model was a key step for the later 
development of dynamical field theory [8]. Another problem in which large d expansion is 
relevant is that of a particle in a random potential. According to the selfconsistent theory 
of localization, [10] explicit analytical results for the critical disorder that induces a metal- 
insulator transition are only known for a Cayley tree geometry which corresponds to a 
lattice of infinite dimensionality. However, there is still qualitative agreement with the 
numerical results in a three dimensional lattice [11]. 

Similarly, many problems in percolation [9] and spin chains [6] have explicit analyti¬ 
cal results in the limit of large spatial dimensions. In many of these cases just keeping the 
leading term in the 1/d expansion is enough to find good agreement with experimental or 
numerical results [3] in d = 3. In the context of quantum gravity, large d expansions have 
also been employed [12, 13] to simplify Feynman diagrams in a spirit similar to the large 
N approximation, broadly used in quantum chromodynamics, N = 3, and other gauge 
theories. However, renormalization of quantum theories of gravity is even more prob¬ 
lematic as dimensionality increases so it is not clear whether it is a viable approximation 
scheme. The situation is different in classical theories of gravity which are finite for any 
dimensionality. The study of properties of black holes [14] and general relativity [15] in 
large dimensions has shown that there are intriguing features that only occur for a suffi¬ 
ciently large number of dimensions. More recently [16-18] this large d limit was studied 
in the context of AdS spaces and then applied, by AdS/CFT techniques, to the study of 
holographic superconductors [16]. 

One of the main conclusions of [16] is that it is possible to find an explicit analytical 
expression of the critical temperature in the limit of large dimensionality and negligible 
backreaction of the scalar on the metric and on the gauge field. Even for d — 2 + 1, this 
simple analytical prediction for the critical temperature is already a good approximation 
of the numerical results. Moreover, as dimensionality increases the condensation of the 
scalar occurs always close to the horizon as the gravitational effects of the black hole are 
only important in this region. 
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In this paper we continue the study of holographic superconductors in the large d 
limit with a twofold motivation. Firstly, we aim to emphasize the usefulness of large 
d expansions in holography by carrying out analytical calculations of the entanglement 
entropy and the conductivity that are only possible in this limit. Secondly, we seek to 
clarify the qualitative effect of dimensionality in holography. We have found that as d 
increases the coherence peak becomes narrower and the ratio between the energy to break 
the condensate and the critical temperature decreases. This is a strong suggestion that the 
effective coupling that controls the interactions of the condensate seems to be weaker as 
dimensionality increases. 

The organization of the paper is as follows: In section 2 we introduced the two mod¬ 
els that we employ to study a large d holographic superconductor, then in section 3 we 
compute numerically the conductivity up to d = 9. Based on these results we compute 
in section 4 the superconducting energy gap, roughly the maximum of the conductivity, 
and the order parameter (O) as a function of d. We also discuss certain ambiguity in the 
relation between these two quantities. In section 5 we study analytically at T = 0, the low 
and large frequency-dependence of the electrical conductivity. Similarly, in section 6, we 
provide simple analytical expressions for the entanglement entropy between a rectangular 
strip and its complement in the boundary; we analyze both the case T = 0 and T ~ T c . 

2 Models 

We study the dimensional dependence of holographic superconductivity [19, 20] in two 
models, one at T = 0 and other at T > 0. 

2.1 d-dimensional holographic superconductivity at T = 0 

For the T = 0 limit we choose the model introduced in Ref.[21], to describe a quantum 
critical point with spontaneous symmetry breaking, 



( 2 . 1 ) 


where 


y(|V>|) = 2A + m 2 |Vf+ |M 4 , 


( 2 . 2 ) 


A = — d(d — l)/2 L 2 is the cosmological constant, nn 2 < 0 is the scalar mass and u > 0. 
Symmetry breaking is directly related to the existence of a minimum of the potential at 


M = V'\R = sj-m 2 /u 7^ 0. 


- 3 - 





Following [21] we consider the metric ansatz 


dr 2 

ds 2 = e 2A(r) (-h(r)dt 2 + d^dx 1 ) + —, (2.3) 

v ' h(r) 

i = 1,..., d — 1, such that in the infrared limit A(r) = r / L\ R and h(r) = 1 where Lir is 
defined through ~ d(d ~ v > = V(|^ir|). 

^IR 

In order to recover the SO(c/-l,l) Lorentz symmetry and SO(d,2) conformal symmetry 
deep in the IR the metric should approach 

ds\ r = e 2r / Lm (— dt 2 + dx l dx l ) + dr 2 (2.4) 


where we have imposed 


h(r) —> h m = 1, A(r) —> ——, as r —> — oo. 

L\r 

Similarly, in the UV limit, the appropriate symmetries are restored provided, 


(2.5) 


h(r) —> /iuv, A(r) —>■ A vw — , as r —* oo, (2.6) 

Ju 

with h\jy and /l uv constants related by the rr component of the Einstein equations: 

{d - \){A!tih + dh 2 A' 2 ) + hV(\ip\) - h 2 ^ 2 - e~ 2A q 2 ip 2 (j) 2 + ^e“ 2 V 2 = 0, (2.7) 


where (j) is the t component of the gauge field. Evaluated at the UV boundary, the previous 
equation, yields 

h\jy = ~j2~- ( 2 - 8 ) 

Moreover, the null energy condition requires huy > h [R = 1 [21] which means that /l uv < 
1. At the same time, A{r) must increase monotonically in the whole range — oo < r < oo 
and the slope in the UV-limit must be lower than in the IR-limit, i.e., 2fuv /L < 1/Lir, 
[ 22 ], 

The resulting equations of motion are, 


r + + -M') + ^ = 0. ^ W - 2)A r - = o, 

..// , ju! AI ^ „2 j.2„/.2„-2A j./2„-2A _ n /t// , ^ „/./2 , e 2A q 2 <f> 2, lp 2 _q 


h" + dh’A ' - -q 2 (j) 2 ^e 
h 


^e~ 2A = 0, A" + 


d — 1 


r c + 


(d~i)h 2 


(2.9) 


with boundary conditions in the IR-limit (r —> — oo), 

$ ~ 0 o e^ [A ^ R_(d_2)1 , ip = + a^e^ { ^- d \ 


( 2 . 10 ) 
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where A^ IR and A, RIR are the larger roots of: A 0|R |A,, IR — (d — 2)] = 2 q 2 ijjf R L 2 R and 
A^, ir (A^, ir — d) = \V"{i^)L 2 lR = — 2m 2 L 2 R . Similarly in the UV limit (r — >■ oo), 

0 = p — pe~^ uvr ; 0 = 0 uv e _ i" 4uv ^ _A ^uv), (2.11) 

where, A dIJV = d — 2 and A, Ruv is the smaller root of: A^, UV (A^ UV — d) = m 2 L 2 /(/r uv Au V ). 
The boundary conditions for h and A are given in eqs. (2.5) and (2.6). Moreover, we will 
take the parameters m 2 and u such that the operators dual to 0 and 0 are irrelevant in 
the IR so that the IR AdS space is a fixed point of the RG flow. Repeating the argument 
presented in [23] it is straightforward to see this corresponds, in our notation, to: 

A^, ir > d, A^ > d - 1. (2.12) 


2.2 d-dimensional holographic superconductivity at T > 0 

For the study of holographic superconductors at finite temperature we employ the, by now, 
standard model introduced in [19, 20] by coupling anti-de Sitter gravity to a Maxwell 
field and a charged scalar and a quadratic (in |0|) potential. Here we state the action and 
equations of motion in d dimensions directly in order to settle down notation and refer to 
the reviews Refs. [24, 25] for more details. The action is given by, 

1 


S = / d d+1 x 


R- F '-\ D ^-Vm 


, R(|0|) = -2A + m 2 |0| 2 , (2.13) 


with T0 = — iff A IL . In probe limit, corresponding to a negligible backreaction of the 

scalar and the Maxwell field on the geometry, is simply given by the planar-Schwarzchild 
AdS black hole, 


r“ . „ , , o L 2 dr 2 


ds 2 = - — h(r)dt 2 + 

L z r z h(r 


+ r 2 dx l dx\ i = 1,..., d — 1, 


(2.14) 


with h(r) = 1 — r^/r d . Assuming for the moment that the only component of the Maxwell 
field is A t = 0(r) it is straightforward to obtain: 

,d -1 20 2 


v ,» +v ,<(* + l±iA v ,A_ + f_M)_° ! 


h" + 


r 2 li 


= 0. (2.15) 


The boundary conditions are fixed from to the usual expansions: 


/ / N a. j3 . p 

^( r ^°°) = 0AT + 0A7 + --- <P(r^oo) = p+-^ + ... 

0(r ->■ r 0 ) = 0o + (l - ..., 0(r —► 7- 0 ) = (l - + ... , 


(2.16) 


where A± = | {d ± V d 2 + 4m 2 L 2 ) and 0i is given in terms of the undetermined con¬ 
stants 00 ) 01 • 
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3 Electrical conductivity in the large d limit: T > 0 case 


We start our analysis by computing the electrical conductivity, cr, at T > 0. For the sake 
of completeness we review the procedure to compute it for a general d. To this end one 
should add a perturbation to the vector potential 5A = A x as well as one to the metric 
$9 = 9tx■ However, we will solve for a numerically in the probe limit where, 

ds 2 — — ( f(z)dt 2 + dz 2 + dx 

z V f( z ) 

with z = 1/r, and choosing the horizon position z 0 — 1 and L — 1. As usual, the 
linear response of an operator, in our case the current J fi (x), to an external source or field 
perturbation, A x , is related, in momentum space, to the retarded Green’s function, [24]: 


f(z) = l-z d , i — 1,... ,d — 1, (3.1) 


fiJ x (k ) = G%{k)5A x (k), 


(3.2) 


where k — (cu, k) is the ('/-momentum and Gff ( k ) is the Fourier transform of the retarded 
Green’s function. Moreover, the charge current response to an electric field is J l {oS) = 

a lj (u)Ej(uj), with E x = —d t A x (t,z,x), A x (t,z,x) = e tut A x (z,x), therefore, it follows 
that, 


<j xx (uj) = 


:«,o) 


IU 


(3.3) 


We now compute this Green’s function following the procedure first outlined in Ref. [26]. 
First, we write the Fourier transform of the vector potential, 


A^( z ,x) 


d d k 

(2ttY 


e ikx A,(z,k) 


(3.4) 


where kx = —cot + k ■ x and A^\k) = A^(z = 0, k) is defined from the boundary value 
A t ,(z = 0, x). The Fourier transform of the gauge-field-part of the action leads to, 


S t 


gauge 


d d k 
( 2 ttY 


■JF(k,z) 


2=1 

+ ..., F(k,z) 


v f=gg zz g x 


A x (z, k)d z A x (z, k), 


2=0 


(3.5) 


where the dots correspond to terms not containing A x and its derivatives. The final expres¬ 
sion for the conductivity is obtained by combining the proposal of Ref. [26] 0) = 

~ 2 6Al°\- S k)sA^(k) lim ^° z ) together with eq. 3.3, 


Refer (a;)] 


1 

iu 5ii°)(-a;)Mi 0) (a;) 


lim ^gg zz g xx A x (z, -k)d z A x (z, k) 


k =0 


(3.6) 
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In order to compute A x (z, k ) we write the equation for A x (z, x ) in the fixed background 
given in eq. (3.1). As was mentioned above, we assume a harmonic time dependence for 
A x . The derivatives are taken with respect to the holographic coordinate, z: 


A'! + 




a: + 





0, / 


l-z a 


(3.7) 


Finally, we impose the usual boundary conditions, in-falling close to the horizon, A x ~ 
(f / z 2 )~ luJ ^ d ( 1 + ...) and A x ~ A^ + A^g d (z, uj) close to the boundary. The function 
g d (z,uj) is easily obtained, for each d, by solving the asymptotic expansion of eq. (3.7), 
App. C. By combining eqs. (3.4) and (3.7) we obtain an analytically solvable differential 
equation for A x (z, k) at zero spatial momentum, which in the z —> 0 limit reduces to, 

Q 

A! x {z, uj) + ^(z, uj) —--b A x (z , uj)uj 2 = 0. (3.8) 

z 

We choose the regular solution at z = 0. We note that at uj — 0 the conductivity Re (A) 
develops a delta function as a consequence of the translational invariance of the system. 

For odd d we have now all the ingredients to compute the conductivity a(uj) (3.6). 
However for even d, logarithmic divergences at non-zero uj appear [27]. In order to study 
the large-<7 limit of a, it is enough to restrict our analysis to odd d. Therefore, in order to 
avoid the intricacies of adding the counterterms to the action to remove the divergences 
mentioned above, we take the prescription for d = 4 given in [27] and for d = 3, 5, 7, 9 we 
employ eq. (3.6). 


3.1 Numerical calculation of the conductivity at low temperature for d < 9 

In this section we compute numerically the electrical conductivity in the probe limit for 
3,4, 5, 7 and 9 dimensions of the dual boundary theory and for two scalar masses m 2 = 
0, d + 1. We follow the procedure described in the previous section and solve the resulting 
differential equations by the shooting method. See Appendix C for the specific expressions 
of the electrical conductivity in each dimension. The results depicted in figure 1 and figure 
2 indicate that as dimensionality increases, the coherence peak becomes narrower and the 
position to the peak u g moves to lower frequencies. The physical interpretation of these 
features is clear. The condensate becomes less coupled as it costs less energy to break 
it (smaller u g ) as d increases. Moreover, the effective bulk coupling also decreases as a 
narrower coherence peak is a signature of a longer life-time of the relevant excitations 
around uj g . A tentative explanation of this behaviour in the gravity dual is that [15, 16] as 
the dimensionality increases the condensation of the scalar gradually occurs closer to the 
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horizon which corresponds to the less strongly interacting limit of the dual field theory. A 
natural question to ask is whether the gravity dual has a well defined limit for d —» oo. 
In order to answer this question in figure 3 we plot oj 9 /T c as a function of d. The ratio 
decreases monotonically as d increases and it is likely to converge to a finite value in the 
d —* oo limit still above the prediction ~ 3.528 of the Bardeen-Cooper-Schrieffer (BCS) 
theory of weakly coupled superconductors. It seems that this limiting value only depends 
weakly on the scalar conformal weight. More specifically we expect this result to hold 
provided that both the chemical potential, related to the kinetic energy, and the conformal 
weight, related to interaction energy, have the same scaling with d. It would be interesting 
to explore whether there exists a strict minimum bound for this quantity in the large d 
limit. 



Figure 1. Conductivity (3.6) in different dimensions for a massless scalar field at T/T c ~ 0.1. 
As the dimensionality increases the coherence peak is narrower and moves to the region of lower 
frequencies. 
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Figure 2. Conductivity (3.6) in different dimensions for m 2 L 2 = d + 1 at T/T c 0.6. Results are 
similar to those of figure 1 for m = 0 

A few comments are in order: a) we omit the case d — 3 for m 2 L 2 = d + 1 in what 
follows since we have observed an anomalous behavior of the AC conductivity similar to 
that reported in Ref.[27], b) in figure 2, corresponding to m 2 L 2 — d + 1, the ‘crossing 
point’ where all curves meet, oj/T c ~ 7, is slightly blurred due to the presence of extra 
poles, not shown, at lower frequency [27], c) although for m 2 L 2 = d+ 1 we found difficult 
to decrease the temperature below T jT c ~ 0.6 it is clear, see figure 3, that the behavior for 
both masses is strikingly similar. 
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d 

Figure 3. Ratio between the peak of the conductivity u g and the critical temperature for m 2 L 2 = 0 
(A = d) at T ~ 0.1T C , and m 2 L 2 = d + 1 (A = d + 1) at T ~ 0.6T C as a function of the 
dimensionality. It always decreases as d increases and only depends weakly on the scalar mass 
m 2 , see figures 1, 2. The horizontal black line indicates the approximate position of the crossing 
point in figures 1 and 2, which we also expect to correspond to the location of the peak of the 
conductivity for d oo limit. We note that even in the d —»• oo limit the ratio is still substantially 
larger than the BCS prediction 3.528. 


4 Relation between the order parameter (O) and uj 9 in the large d 
limit 

In BCS superconductors the coherence peak in the conductivity is simply two times the 
value of the order parameter also referred to as the superconducting energy gap. Physically 
it means that since a Cooper pair is composed of two electrons it takes twice the energy 
gap to break a Cooper pair and place these two electrons in the first state available above 
the Fermi energy. For strongly coupled superconductors there is no clear relation between 
these two quantities as the coherence peak broadens substantially and in some materials 
the quasiparticle picture based on the Fermi liquid approximation breaks down. However, 
in the context of holographic superconductivity it is well known [28] that these two observ¬ 
ables are still comparable, though the relation between them is not universal and different 
from the BCS prediction [29]. We now study to what extent this relation still holds in the 
large d limit. The order parameter ( O) is computed by following the usual steps. First we 
find the numerical solution of the equations of motion eq. (2.15) by the shooting method 
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for a scalar field, ip(r), charged under the gauge field A, = Mr) in a non-dynamical 
Schwarzchild background, i.e., in the probe limit. We consider a fixed charge density and 
different scalar masses. The order parameter is simply ( O ) A = [(2A — d)3\ A , where A 
is the conformal dimension of the operator dual to t/>, d is the number of dimensions of the 
dual theory, and [3 is given in the boundary condition, eq. (2.16). In Table 1 we present 
results for ui g and (O') for different dimensions and masses. As dimensionality increases 
( O) becomes much smaller than u g . Indeed, it seems that the ratio (0''M A /oj g —> 0 as 
d —> oo. Presently we do not have a solid explanation for this discrepancy. A finite value 
of the order parameter ( O) in holographic superconductivity is interpreted as a signature 
of spontaneously symmetry breaking rather than a energy gap in the spectrum. It might 
therefore be that these two quantities are not related and the similar value in low dimen¬ 
sions is a coincidence. Another more speculative explanation is that the standard recipe to 
compute ( O) misses some dimensionality prefactor. We went over the original derivation 
of the expression for the order parameter but we could not find any discrepancy with the 
expression used above. However we found that by rescaling, see figure 4, (O) by T(A) 
the ratio seems to converge to a finite positive value in the d —> oo limit. 1 Whether this 
is just a coincidence or has a deeper physical meaning remains to be understood. Finally, 
we note the fact that the rescaling by T(A) depends on the scalar mass indicates that it is 
not related to the dimensional dependence of the coupling constant in the action which is 
usually set to the unity. 



UJg 

T c 

( O )V A 
T c 

(0)!/A 

T c 

(o> 1 / A 

UJ q 

<i 6 > J / A 

UJq 

d = 3 

11.3 

12.7 

16.0 

1.1 

1.4 

d = 4 

10.0 

8.7 

13.6 

0.9 

1.4 

d — 5 

9.1 

6.4 

12.1 

0.7 

1.3 

d = 7 

8.4 

4.1 

10.5 

0.5 

1.3 

d = 9 

8.0 

2.9 

9.4 

0.4 

1.2 


Table 1 . Comparison of the position of the conductivity (3.6) coherence peak with the order pa¬ 
rameter (O) = (2A — d)/3 and the alternate definition (O) = (2A — d)T(A)/3, for m 2 L 2 = 0 
(A = d). Convergence for large d is only observed after the order parameter is rescaled by T(A). 
We do not have a clear understanding of why the order parameter and u g have a different parametric 
dependence on the dimensionality. 


'Our numerical results suggest [(2A d)rf\ a —> constant for d oo. Thus, a factor depending only on 

d such as T(d), instead of A, does not result in a finite (O) 1 ^/T c in the the limit d —> oo. 
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Figure 4. Ratio between the order parameter close to zero temperature and the critical temperature 
for different dimensions and scalar masses. Dashed lines correspond to the usual definition: (G) = 
(2A — d)/3, continuous lines include a speculative factor T(A), (O) = (2A — d)T(A)/3. Only in 
the latter case convergence of the ratio to a non-zero value in the d —> oo limit is likely. 


5 Analytical calculation of the conductivity at T = 0 for different di¬ 
mensions 


We now switch to the background introduced previously in section2.1 to describe holo¬ 
graphic superconductivity at T = 0. From now the main focus of the paper will be to 
compute analylitically the conductivity and later the entanglement entropy in the large d 
limit in order to illustrate the interest of large 1/d expansion in holography. 

In this section we compute the electrical conductivity at zero temperature. As was 
mentioned previously one must consider fluctuations of A x (t, r ) and g tx (t, r ), [20], which 
source an electric field E x and carries momentum T tx . These perturbations are usually 
assumed to have a harmonic time dependence, A x {r)e ~ luJt , gi.x(r)e^ VjJI '. Furthermore, the 
Einstein and Maxwell equations are expanded in g tx (r) keeping only linear terms in A x {r), 


A"(r) + A'(r) 


{d- 2)A' + ^ 
h 


+ 


A x (r) 

h 


U! 


he 2A 


2gV - (j) ,2 e~ 2A 


= 0. (5.1) 
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We then impose that near the UV boundary A x (r) = A 0 + T > A(r >. In the infra-red 

limit we expect the perturbation A x to become small. This is indeed the case for d = 3 
but not for d > 4 [30] where it grows exponentially for r —>■ — oo. This cast doubts about 
the stability of the background to small perturbations in large dimensions. Indeed, it has 
been observed that the addition of a gauge field increases the temperature of the dual field 
theory [31] even in the limit of an extremal black hole. However a full stability analysis is 
beyond the scope of the paper as the main motivation here is to employ the large d limit 
as a computation tool to obtain analytical results. As in the d = 4 case studied in Ref. 
[32] we overlook the potential instability induced by the gauge field and proceed to solve 
analytically eq. 5.1 in the following three different limits. 


5.1 Low frequencies 

The small frequency dependence of a is studied by solving eq. (5.1) in the IR limit. The 
scalar is now locked around its minimum, if i R . By using the asymptotic values of A and h 
in the IR limit eq.5.1 simplifies to, 

A xi r ) + + A A r ) ^ = o, (5.2) 

2r 

where we have assumed that e 1 ® <j)'(r) —» 0 as r —» — oo. The solution of the above 
equation can be written in terms of a Hankel function as: 

A x {r) = (uL lR e~A^j , a = A^ m - ^—1 = I ^(d - 2) 2 + 8q 2 ipf R L 2 R . 

(5.3) 


As was pointed out previously, [21, 28], the frequency dependence of the conductivity 
at zero temperature is extracted from the conservation of the flux (d r T = 0) with T = 

—fte (d - 2)A A 

2i J • rL X9 


R e (a) oc ^ (5.4) 

u\A 0 \ 2 

Notice that, modulo a factor if 2, the flux T coincides with the definition of J-(k, z ) 
given in eq. (3.5), namely ^(k, z ) = — 9 A x {z, —k)d z A x (z, k ), where in this case 

the metric is given by eq. (2.3). In the latter the holographic coordinate is r, instead of 
z — 1/r, and we take the gauge field in position space instead of momentum space. 

To obtain A 0 we need to match the solution given in eq. (5.3) to Z(r), the solution of 

_ r_ / cl—2 \ 

eq. (5.1) with u = 0, which is assumed to satisfy Z(r) —> e 2 > as r —>■ — oo. 
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For u small enough, such that r* <C r IR , where r* = L IR logo;L IR and r IR is the scale 
at which the geometry is significantly deformed from eq. (2.4), the convergent part of the 
solution, eq. (5.3), is matched to Z(r) in the region r* -C r -C ri R : 

A x (r)~C Z(r)(uL m )- a , (5.5) 

where C is a constant. Therefore, taking the limit of the previous expression when r —* oc 
results in 

A 0 oc co~ a , (5.6) 

and, from eq. (5.4), the conductivity is, 

Re (a) oca; 2 "” 1 , (5.7) 

with a = A 0UI — ^=2 — \yj{d — 2) 2 + 8 q 2 ^ R Lf R . The exponent a that controls the 
strength of the low energy excitations increases with d. This is a strong suggestion that, 
in agreement to the results at finite temperature, high dimensionality suppress low energy 
excitations and therefore make the system less strongly interacting. The d dependence of 
the conductivity in low frequency limit was previously investigated in Ref. [30]. However 
the expression for the conductivity in [30] is not the same as eq.(5.7). We note that eq. (5.7) 
agrees with the results of Refs. [21, 32] for d — 3,4 as well as with our numerical results up 
to d = 9. We observe that as d increases, the region where eq.(5.7) is a good approximation 
is restricted to smaller frequencies. Moreover, since for larger d the divergence in A x is 
stronger, see eq. (5.3), the numerical results become less reliable, and harder to obtain, in 
this limit. 

Finally we also note that, in Lifshitz backgrounds with hyperscaling violation, the DC 
conductivity for small frequencies shows a similar power law behavior [33, 34], It would 
be interesting to carry out a 1 /d expansion in these type of backgrounds in order to explore 
universal features in the large d limit. 


5.2 Large frequencies 


We now explore the large frequency limit of the conductivity corresponding to the region 
where the frequency u is the largest energy scale in the problem, namely, it is much larger 
than the chemical potential or the condensate. Since the conductivity has units of energy 
to d — 3 we expect that in this limit its real part oc cu d-3 . This can be confirmed explicitly 
by rewriting the prefactor in front of A x in the third term of the left hand side of eq. (5.1) 
as: 


l (r e ~ 2A ~ 2g2 ^ 2 _ 2g " 2A ) • 


(5.8) 
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For ujL —> oo and r —>■ oo such that uLe r / L ~ 0(1), the last two terms are negligible 
with respect to the first, 


2q' 2 J(r) 2 + e 2A ^J(r) 


^Q 2 ' l l J uv' 


-~z A uv ( rf - A % y ) _|_ 


p 2 (d 


L 2 


2) 2 

-e 


2d 1 


(5.9) 


by virtue of the boundary conditions given in eq. (2.11). These terms are negligible 
compared to the term OC LU 2 , which by assumption is 

~ 0(1). (5.10) 

in the region of r considered. For even larger r, the previous term becomes arbitrarily 
small and no additional u dependence is introduced. Thus, all the frequency dependence 
of A x is obtained, in this region of frequency, by setting the scalar to zero and h ~ 1 and 
A ~ r/L. This leads to 


, , , d-2 

A Jr) = e 


H% (coLe~t) + CJH% (uLe~Z) 


( 2 ) 


C*2 G 


(5.11) 


C ‘2 has to be determined from the solution in the bulk, however since we set h ~ 1 and 
A ~ r/L in the whole domain of r, the solution in the IR is approximatively given by 
setting the solution above. Therefore, the ingoing boundary conditions imply C 2 ~ 0. 
Physically, for large enough u the perturbation is insensitive to the flow between the two 
AdS spaces, and, in particular, to the presence of a nonzero scalar field in the degenerate 
horizon. 

Close to the boundary, eq. (5.11) reads 


A Jr) 


2-d 

Clu 2 


+ ... 


(5.12) 


where C is a constant and the dots stand for terms which depend on uj but decay exponen¬ 
tially with r. Therefore, 

A 0 = lim AJr) oc (5.13) 

r—>00 

and, as before, using eq. (5.4) leads to, 

Re (a) ~ CJ l u d - 3 (5.14) 


with 

C d = n -[T(2-d/2)]- 2 2 d - 2 . (5.15) 

We note that this result is strictly valid for odd dimensions only. The case d = 4 has been 
discussed in [27] where, it was found a ~ u[n/2 + ij + log ^-)] for large frequencies. 
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6 Analytical calculation of the entanglement entropy in d /$> 1 di¬ 
mensions 

The entanglement entropy is a valuable source of information of strongly interacting sys¬ 
tems including the classification of the different quantum phases, the estimation of the 
effective number of degrees of freedom of the theory, the rate of propagation of informa¬ 
tion after a perturbation or the location and characterization of phase transitions even in 
cases where there is no order parameter. In the context of holography it has also been 
intensively investigated after the landmark conjecture of Ref. [35] provided a relatively 
straightforward procedure to compute it. Several papers [36-40] have already discussed 
the entanglement entropy in holographic superconductors [37, 38], metal-superconductor 
transitions [39], metal-insulators transitions [36] or in a superconducting interface [40]. 
It has been found that the entanglement is a good observable to characterize these transi¬ 
tions. Its value is always smaller in the condensed phase and has a discontinuity or a kink 
(discontinuous derivative) that signals the transition point. It is also sensitive to a mass 
gap or to the proximity effect in an interface. These calculations in holographic supercon¬ 
ductors are numerical as the calculation of the entanglement entropy requires to compute 
the backreaction of the scalar and gauge fields on the background. The main goal of this 
section is to show that explicit analytical results are possible in certain T = 0 backgrounds 
and also around the critical temperature but only in the limit of large spatial dimensions. 
This is a strong indication that 1/d expansions in holography broadens substantially the 
scope of the problems that can be addressed analytically. 

6.1 Entanglement entropy at zero temperature 

We now calculate analytically the entanglement entropy at zero temperature related to 
the background eqs. (2.1)-(2.6). According to the usual prescription [35] proposed by 
Takayanagi and Ryu, given a field theory in d dimensions, the entanglement entropy of a 
region of space A and its complement is calculated from the gravity dual by finding the 
minimal d — 1 -dimensional surface 7 ^ which extends into the bulk such that O 77 = <9/1. 
In other words, the boundary of 77 at the AdS^+i boundary is equal to the boundary of A. 

To illustrate the calculation we choose A to be a d — 1 dimensional strip of width t. 
I = {i6 M d_1 : — 1/2 < x 1 < 1/2, —a/2 < x l < a/2, i = 2, ..., d — 2}, where a is 
the “length” of the strip. The entanglement entropy related to the metric eq. (2.3) is, 



( 6 . 1 ) 
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where, x = x 1 , a d 2 results from integrating x l with i = 2 ,..., d — 2 and Cff 1 is the 
d + 1-dimensional Newton’s constant. 


6.1.1 Sharp domain wail approximation 

As was mentioned above, the background given in eq.(2.3) interpolates between two copies 
of AdS space in the IR and UV regions, eqs. (2.5) and (2.6). Since there is no analytical 
expression for h(r) and A{r) in the whole range of r we follow [41] and assume a sharp 
transition between the two AdS domains at a position denoted by r DW - Numerical results 
show that there exists a — oo < r m < 0 such that 'w'{r rn ) =0. It is therefore natural to 
choose r DW = r m . Even though we will not be interested in the specific value of r DW 
we will require r DW < 0 in the following sections. Moreover, numerical results suggest 
r m oc dr 1 . 

More specifically the sharp domain wall approximation consists in taking A{r) and 
h(r) as the asymptotic values given in eq. (2.5) for r < r DW . Similarly, for r > r DW we 
take those given in eq. ( 2 . 6 ). 

As usual in the calculation of S 7 with A a strip, eq. (6.1) does not depend on the 
integration variable x explicitly. Therefore, the Euler-Lagrange equations that minimize 
reduce to the Beltrami identity which states that, given a Lagrangian L, if dL/dx = 0, 
then L — r'dL/dr' is a constant. In our case: 


e {d-l)A{r) _ J e (d-l)A vw ^ , r> T DW 

y/l + e~ 2A ^r' (x) 2 / h{r) y e^ _1 ' >x iR ; r < row 


( 6 . 2 ) 


In the previous equation we took into account the different AdS radii, L and L IR in each 
region, and r* is the “turning” point of the surface 7 ^ which occurs for x = 0. We will 
consider the general case r* < r^w, he. the minimal surface extends into the IR region. 
With the previous considerations eq. (6.2) is easily integrated, 

£ 

dx = - = Jir + Juv, (6.3) 

where 



rr dw 


hlR — 


dr 


m jv 


e^iR 


hsv — 


dr 


2(d-l) 7 - 


- 1 


'nsw 


v//q^e Auv i x / e 2(d-1)A uv 


r — r* 
L 


(6.4) 
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r uv being the UV cutoff, Jir is calculated with the change of variables t — e 


2(d—2) - 


Iir — ~Li R e 




d nm) 


nr + 


2 F 1 


d 3d — 2 


2(d—1) - 


2’ 2d — 2’ 2d — 2’ 


(6.5) 


while an analogous change of variables, t = e 2(d 2,/l ' JV S w in Juv yields 


Juv — L- 


; [(d-l)r«-dr D w]- 

d 


2^1 ", 


d 3d — 2 


2 2d — 2 2d — 2 


,e 


2(d—l)Auv 


r * — r DW 


( 6 . 6 ) 


where we used the relation between d uv and /l i;v given in eq. (2.8). Similarly, inserting 
eq. ( 6 . 2 ) into eq. ( 6 . 1 ), the entanglement entropy can be calculated by integrating in r in 
the two domains (r < r DW and r > r DW ): 




2 a d ~ 2 

AG d+ 1 [*^ IR + ^uv] , 


(6.7) 


crow 


Sir — 


(d—2) 7 

r l ir dr 


T (d-2)f*. 

L i R e 


—2(d— 1 ) rzi* 2(d — 1) 

6 ^IR 


dw 


y 


1/2 


w 


1/2 


(1 — yu) 2d ~ 2 


T (d-2)f±- 

SiRe ;l ir 
2(d — 1) 


( 6 . 8 ) 


n ^ n . 1 3d-4 3 
1 ^ 2’ 2d - 2’ 2 ,2/ 


where we made the change of variables: ti(r) = 2 (l — e - 2 ( d - 1 )(»’-» i .)/ i iR) ) y — 1 _ 
e - 2 (d-i)(r' D w-r*)/L IR C q (g gj ca n |-) C rewritten using the following Hypergeometric func¬ 
tion identities, 


2 Fi(a, b , c, z ) = r ^ ~ T~Ty\ 2-Si( c - a, c - 6, c - a - b + 1,1 - z) + 

1 (1 — c)i (c — a — b + 1) 

r(i — a)r(i — 6)r(c) . 

1 -z 2 Fi(a — c + 1 , b — c + 1 , 2 — c, z), 


r(2-c)r(c-ia)r(c-6) 


(6.9) 


with 6 = 0 , c = 1^1, a = c — 1/2 and 

c 2 Fi(a — 1, 6, c, 0 ) — c 2 Fi(a.6 — 1, c, z) — (a — b) 2 i 7 i(a, 6, c + 1, 0 ) = 0, (6.10) 


with a — 1 / 2 , b = c = as follows 


Sir —Sir 


(d—2) jz- p/ 3d^2 \ (d_2)32tt 




ilRi( 2d=iJ , e ^.r 1 / i _ e _ 2(d _i } !m=n + 


d(d — 2) r(|5i) + d — 2 
r ^+ 2 (d-r)^ / X d 3 d -2 


d(d — 2) 


2^1 ^ 


2’ 2d- 2’ 2d- 2’ 


r»-r DW 

e ' ' ^ir 


^2(d—1) - 
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On the other hand, for ,5' uv , one must take care of the usual divergence for r —> oo. 
Defining the auxiliary variables t = e~ 2 ( d_1 ) Auv i DW , the cutoff in the t variable f uv = 

e -2(d-l)A U v DV i D * anc l y — e - 2 (d-l)A vwwe j nte g ra( - e ^ uv; 


S'uv — 


l nr uv e (d-2)A uv £ d r 

y/hirv J r DW y/l — e -2A uv (d-i) 

Le (d-2)A m r AM- 


Le (d-2 ) Avy nm ,1 dt 


i —r* 

L 


2 (d - 1) J tvw JmJ y/(l - yt) 


= L 


2{d — 1 ) 

g(d— 2)A v -v^X g(d—2 )A uv ^j^ 


q--= rr - /I 3d — 4 3 

22 / 2(d_1) \Jl — ty 2 Fi —— 3 , 3,1 -yt 


d-2 


d-2 


3 [-dr DW +2(d-l)r*]^ff 

d(d — 2) 


2*1 n, 


2’ 2d — 2 2 

^ _ g—2(d-l)Auv rp \ r * _[_ 

d _ 3d - 2 g2 (d-i)A uv T '*~ 6 rDy 


t=i 


- t=tuv 


2’ 2d- 2’ 2d- 2’ 


( 6 . 12 ) 


In the last equality we used the relations given in eqs. (6.9) and (6.10) and left the cutoff 
r uv explicit in the divergent term. 

We stress eqs. (6.7), (6.11) and (6.12) are an approximation to the entanglement 
entropy between a strip of width £ and its complement in the d-dimensional boundary when 
the scalar field condensates. It is interesting to compare these results with the entanglement 
entropy between a strip of the same width £ and its complement in the situation in which 
the scalar is absent, [42]. To do so we should express S^ in terms of £, however, from eqs. 
(6.3), (6.5) and ( 6 . 6 ) it is clear that r* cannot be expressed in terms of £ in a closed form 
and thus the comparison cannot be made easily. Instead, in the next section we make this 
comparison only in UV and IR limits of S 3 . Additionally, we also study the large-d limit 
of% 


6.1.2 UV, IR and large-d limits 


UV limit: we first consider r* > r DW , he. the minimal surface 73 is embedded in the AdS 
copy that contains the boundary r —>■ 00 . In this situation Ji R = 5 )r = 0 and r DW = r* in 
eqs. ( 6 . 6 ) and ( 6 . 12 ): 


cuv 

^A 


rv_/ 



\ = Le ~ 
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As was expected we recover the result for the infinite strip in an AdS space, found in 
[43]. It is observed the strip width tends to zero following e~ / ' uv v\ while the “finite” part 
of the entanglement entropy diverges as e ( d - 2 h 4 uvI r. 


IR limit: In case r* -C r DW , i.e. the 77 extends deeply into the IR region. From eqs. 
(6.5), ( 6 . 6 ), (6.11) and (6.12) 
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Li R e L[R 


v^r(jEi) 

d r mv 


oIR 

^4 


rsj 


2 a d ~ 2 

'^ e (d-2)A v Gf /2\ d - 2 L?- 1 ^ 
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d- 2 \i) d- 2 

[ r ( M - 2 )J 

, Li R e( 

- 2 )^“ _ Le (d-2)A VW ^ ] 
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d- 2 j ‘ 



d— 1 


+ 


(6.14) 


In this limit, the strip width, l, diverges and the finite part of the entanglement entropy 
saturates to a constant value given by the first term in the following expression: 


lim Sf 


2 a d ~ 2 

4G d N + \d-2) 


L lR e (d ~ 2) 


r DW 

L IR 


Le (d-2)A uv ^w 


2 a d ~2 

+ 4 G^ 1 ~d--2 


(6.15) 

At this point, it is easy to compare eq. (6.14) to the entanglement entropy between the 
strip A (of same width £) and its complement in the d-dimensional AdS boundary when 
•0 = 0 . Were the scalar field be zero, there would be a single AdS space and would 
be given by the first two terms of eq. (6.14) while the last term would be zero for all 
d. In the presence of the condensate, the third term in the previous equation corresponds 
to the contribution of the flow from one AdS copy to the other. Indeed it is easy to see 
that this term is negative. From the definition of Li R : —d(d — l)/L‘f R = R(|0i R |) = 
— d(d — 1 )/L 2 — m 4 /(2w), it follows L IR < L. Since we require r DW < 0, the term in 
square brackets of eq. (6.15) 


L K e‘ d - 2 » 


r DW 

L IR 


Le ( d -2)A UV ^ < e(rf -2M uv r® (Lir _ L) < 0 . 


(6.16) 


The conclusion is that the entanglement entropy between a strip of length t and its 
complement is lower if the scalar is present. This means the theory has less degrees of 
freedom in this case. In the limit of a strip of infinite width (£ —> 00 ), the finite contribution 
of S~ x reaches the maximum value given by the first term in eq. (6.15). These results are 
expected as the entanglement entropy counts the degree of freedom of the theory. It is 
therefore natural that it is smaller in the condensed phase. 
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Let us turn to the study of the large-d limit of S^. Before we do so, we must analyze 
the behavior of the strip length as d —* oo. From eqs. (6.3), (6.5) and (6.6) it is clear that if 
r* either remains constant or increases, as d increases, both Ji R and Juv would tend to zero 
and / —» 0. In order to compare S for different d we must keep £ constant. Therefore, 
r* —>- —oo, as d —> oo, which corresponds to the IR limit (r* <C r DW ) studied above. 
Taking d large and £ constant in eq. (6.14) yields, 


and 


r*(d —>■ oo) 
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(6.18) 


The second term of the previous equation corresponds to the universal contribution for the 
infinite strip in an AdS space, [43] which is strongly suppressed in the d —>■ oo limit as 
it is proportional to d~ d . The third term has some interesting features. It is independent 
of i as it is expected in gapped systems where the typical length, the numerator in this 
case, is closely related to the coherence length of the holographic superconductor. Its 
('/-dependence, arising from the flow of one AdS space to another, is dictated by the in¬ 
dependence of r DW and A uv . As mentioned before, numerical results for d < 9 suggest 
oc d~ l and r DW < 0 which implies a behavior like d~ l for these contributions. In the 
limit of a vanishing scalar field, the third term vanishes for all d and the result of Ref. [43] 
is recovered. 

Let us simplify eq. (6.18) for the particular set of parameters: m 2 L 2 = —3c/ 2 /16, 

, L — 1. These values, together with the definition of L m : 
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m 2 L 3 
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= V(IV’irI), yield a constant, in d, Li R = ^/32/35 L. Moreover, as discussed 


earlier, -p < -^L. These considerations allow a further simplification of eq. (6.18), 
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(6.19) 


where i and a are the width and the characteristic length (infinite) of the transverse dimen¬ 
sions of the strip. The radius of curvature of the IR asymptotic AdS space, L\r. does not 
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depend on d, a is the constant of proportionality in r DW ~ which can only be obtained 
numerically and A L = L IR — L= \j32/33 — 1 < 0. S dlv is the (divergent) part containing 
the UV integration cutoff. 

As we mentioned previously, were the condensate vanish, the last term in eq. (6.19) 
would be identically zero, since the asymptotic IR and UV AdS radii would be the same, 
L m = L. Moreover, this term is negative, which means the finite part of the entangle¬ 
ment entropy is lower, and thus indicates less degrees of freedom in the presence of the 
condensate. Finally, as d —* oo, this contribution is smaller, suggesting the difference be¬ 
tween the entanglement entropy in the presence and absence of the condensate is smaller. 
The latter is an indication that, in agreement with the conductivity results, the condensate 
interactions become weaker as d increases. 

6.2 Entanglement entropy close to the transition 

In this section we compute analytically the entanglement between the semi-infinite strip, 
A, defined in the previous section and its complement at finite temperature. We employ 
the action eq.(2.13) but we have to go beyond the probe limit. We assume the following 
parametrization of the metric: 

T 4 

z)e~ x ^dt 2 + ——— dz 2 + dx 

/(*) 

where i — 1,... ,d — 1, z — 1/r. Above the transition, the metric corresponds to the AdS 
planar Reissner-Nordstrom in d + 1 dimensions. More precisely, \{ z ) — 6- the gauge field 
A t = <t> = n[ 1 - ( z/z 0 ) d ~ 2 ] and f(z) = /rn = 1 — (1 + Q 2 ) (^J + Q 2 (, where 
Q 2 = // 2 Zq 7 2 , y 2 = and z 0 is the inverse of the outer horizon r 0 . 

Throughout this section we take d to be large so we can get explicit analytical results. 
We also consider strips of length t for which the minimal surface 77 , associated to the 
strip, does not extend too deeply into the bulk, such that the turning point, 27 , satisfies 
{z^/zqY <A 1. This is in general a good approximation in the d —» 00 limit, even for 
z* < Zq. Moreover we restrict ourselves to the region T ~ T c and therefore, the dual 
order parameter (O) is very small compared to the typical energy scale T c . This regime 
restricts the generality of the results for the entanglement entropy but allows to estimate 
analytically the correction in the presence of the scalar field close to the phase transition. 


( 6 . 20 ) 
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The entanglement between the strip and its complement is given by: 
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where we have rescaled z — > z/L 2 in order to compare with the results in Ref. [37]. We 
have also introduced the UV cutoff z\jy —> 0 and, as before, we have used the fact that 
the integral does not depend on x. The turning point, z*, of the surface 7 ^ embedded into 
the bulk is given by z’^ 1 = z d ~ l \/l + z'(x) 2 / f(z). The strip width, £ is related to z* as 
follows: 


„d—1 


dz 


( 6 . 22 ) 


'f(z) (z 2d ~ 2 - z 2d ~ 2 ) 

Even in the absence of the scalar field in eq.(2.1), i.e., the Reissner-Nordstrom back¬ 
ground, the previous two integrals cannot be computed analytically for arbitrary d. How¬ 
ever, an analytical calculation is possible in the large d limit. 

First, we calculate the width of the strip from eq. (6.22), by setting f(z) = Jr\(z) 
and expanding y/ f-m(z) in powers of z/z 0 : 2 

0 „ r ( 3d ~ 2 ) oo n /„\Kl 
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2 n (n-l)\l\ T(^) a nt + 1 ’ 
a = 1 + Q 2 , (3 = Q 2 and a nl = 2 dn + d(l — l) + 2(1 — n) — 1 , b n i = (2d — 2)(n — l)+ dl. 

In the large d limit, assuming t fixed, it is enough to keep only the terms corresponding to 
n — 1 in the series above. The resulting expression of the strip length, l, as a function of 
the turning point of the minimal surface, z*, is, 
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Similarly, from eq. (6.21), with f(z) = / RN (z), 
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2 For simplicity it is more convenient to expand in 6 = — (1 + Q 2 ) (j^Jj + Q 2 
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where C ni is given in eq. (6.24), a n i = (2d—2)(n—l)+d(l—l)+l,b n i = (2d—2)(n—l)+dl. 
For large d, taking the first two terms of the series, 


^ 2 _ 7r 1 + Q 2 fz*\ d _ nQ 2 f z *\ 2d 2 

dz^ 2 d?zt~ 2 2zt~ 2 V^o/ 2zi~ 2 V^o/ 

~ 2 T&-' 1 1 + Q*d? e Q 2 dV 

dzy Y 2 d d i d ~ 2 2 zq 7r 2 27r d - 1 Zo d 2 


where, in the last equality we substituted z* ~ which is a good approximation as 
long as z* -C zq (small £) and d fixed or, for a fixed z* < z 0 , and sufficiently large 
d. In the latter case, £ should also be small, which means that as d increases the minimal 
surface A should reach the near-horizon region for smaller strip lengths. 0 is related to the 

chemical potential and the position of the outer horizon, through // and z 0 , Q 2 = [rz^/ 2 , 

2 _ d^2Lt 
I d—1 2 ' 

In the presence of the scalar field, -0, analytical results are harder to obtain close to 
the phase transition T <T C since f(z) is subject to the backreaction of U', and, in general, 
cannot be written in a closed form. 

However, we show below that it is still possible to find an explicit analytical expression 
in the large-d limit. 

In order to proceed we solve perturbatively the equations of motion close to the 
transition. To do so we expand the fields in the equations of motion (see the appendix 
A for more details) in a power series in a quantity related to the VEV of the operator 
dual to the scalar field. More specifically, from the UV boundary condition of the scalar 
field, -0 ~ + ..., given in eq. (2.16), we set a = 0 and define e = fi. 

Close to the transition this expansion parameter is related to temperature in the usual way, 
e A+ oc (O) oc (T — T c ) 1//2 , with A + being the conformal dimension. 

The blackening function can be expanded as f(z) ~ /rn + e 2 (/f (z) + ...) with f%(z) = 

‘ / \ d / \ 2d-2l 


— HokZq yzo) ~ j an 6 the dots indicate subleading terms, see appendix eq. 
(A. 10), where // ( , is the chemical potential at the phase transition and k is an integration 
constant which is calculated from the perturbative analysis of the equations of motion, eq. 
(A. 12). It is negative k < 0 for d > 3. 

The calculation of the entanglement entropy including the leading correction e 2 /f (z) 
is totally analogous to the one corresponding to the Reisnner-Nordstrom case given in 
detail above. The main difference is that a and /3 in eqs. (6.23) and (6.26) are replaced by, 


i , r\ 2 i 2 ~2 o r\2 , 2 ~2 

Oi — 1 + Q +6 /J,qKjZq, (j — Q + t /IqIZZq. 


(6.28) 












Here, Q 2 = HoZqI 2 ^ Q 2 and 2 0 ^ z 0 , in order to take into account the different horizon 
radius with respect to a pure Reissner-Nordstrom black hole at the same temperature. 

Consequently, in the large-d limit, the relation between the strip width and the turning 
point of 77 in the hairy black hole background is, 


7T a I z* 

2 + 2 d\Y 0 


/3tt 

8 d 


2d—2 


Similarly, §a is 




7T 1 + Q 2 + fl 0 Ke 2 ZQ 


rl ~ d ^ 

U//6^jy 


d 2 zl 


d—2 
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22: 


d—2 


— I - 7T 

Z 0 


+ ... 


Q 2 + n.QKt Zq f 2 


(6.29) 


22 : 


d—2 


Z 0 


2d—2 


+ ■ 


(6.30) 


where, k < 0 is given in the appendix A. In order to compare the entanglement entropy 
between the strip and its complement in the condensed phase with the one in the symmetry 
unbroken phase one needs, in principle, to compute the charge, Q, and horizon position, 
z 0 , of a Reisnner-Nordstrom black hole at the same temperature, eq. (6.31). However it 
is important to note that the contribution due to the condensate, contained in the e 2 term, 
always leads to less entanglement in the condensed phase (/; 0 > 0 and n < 0). 

To compute the Reissner-Nordstrom black hole parameters at the same temperature 
as the hairy black hole we fix the horizon in the superconducting phase, z 0 = 1, and solve 
the following equations in the horizon, z 0 , and charge, 0: 

Q T 1 d-{d-2)Q 2 

d -i i , i i (6.31) 

Zo 7 T 4tt (Q/ 7 )3=t 

where T/p^ i s a function of e (proportional to (0) l ' A ), the metric components at the 
horizon and the chemical potential at the the transition, // 0 , eq. (A. 13), and // = // 0 + 
€ 2 (k + (j) 2(0)) > Ho, where k and 0^(0) are integration constants given in the App. A. 
From eq. (6.31) it follows that z 0 > z 0 and Q < Q, and solving eqs. (6.25) and (6.29) 
one obtains 2* > z*. Consequently, comparing the finite contributions in eqs. (6.27) and 
(6.30), we conclude that below, but close, to the phase transition the number of degrees of 
freedom in the dual field theory is smaller than in the normal phase (no condensate, e = 0, 
/i = //()). This is again consistent with the theoretical expectation that the entanglement 
entropy is closely related to the effective number of degrees of freedom of the system at a 
given temperature. 

For completeness, we express s ; x in terms of the strip length i, the expansion parame¬ 
ter e oc (0) 1/,A , po and k. From eq. (6.29), 2 * ~ ^ [l + 0(d~ 2d )]. Substituting z* = ^ 
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in eq. (6.30), the final expression of the entanglement entropy in terms of the strip length 
is given by, 



cf r +.... 


(6.32) 


Before we compare eq.(6.32) with numerical results we discuss the limits of applica¬ 
bility of the linear approximation z* ~ 2* oc l. 

In figures 5,6 we depict the dependence of the tip of the surface A on £ resulting from 
solving eqs.(6.22) and (6.23). For small i, the linear approximation agrees well with the 
exact result, eq. (6.22). As £ grows this agreement worsens substantially. Additionally, 



is valid for smaller values of t 


but, at the same time, since the corrections 0(d ~ 2d ) are smaller, it remains a good approx¬ 
imation closer and closer to the horizon for both the normal (T < T c ), figure 5, and the 
condensed phase close to the transition, figure 6. This is nothing else but a consequence 
of the simplification of general relativity in the large-d limit. For a black hole, as dimen¬ 
sionality increases, its region of influence shrinks to a neighbourhood of the horizon [15]. 
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t/L 


Figure 5. Position of the tip, 2*, of the minimal surface in the Reissner-Nordstrom background for 
d = 3 and d = 6 . For d = 3, /r = 2.02, Q 2 = /i 2 2o7 2 , 7 s = 331 and z o = 0.992 (dotted line), 
while for d = 6, /1 = 0.38 and 20 = 0.988 (dashed line). The ’’numerical” results are obtained 
from the numerical integration of eq. (6.22) with f(z ) = /rn- The linear approximation 2* oc £ 
corresponds to n = 0 in eq. (6.23) and a = 1 + Q 2 , 3 = Q 2 . The analytical solution of the 
fourth degree polynomial in 2* contains the leading correction, the first term of the series given in 
eq. (6.23). The linear approximation 2* oc i is clearly only valid for small i and, for larger d, it 
becomes gradually better deep in the bulk. Including more corrections in higher powers of 2 * /20 
gives a better approximation but requires, in general, numerical methods. 
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Figure 6. Position of the tip, z*, of the minimal surface in the superconducting phase for d = 3, 
/x = 2.00 and d = 6, /i = 0.37. In both cases the horizon is fixed at zq = 1 (dashed line) and 
Q 2 = /Xq5o7 2 , 7 2 = tt- Similarly to figure 5 the numerical results are obtained from eq. 
(6.22) with f(z) = /rn + e 2 /2 and fi given in eq. (A. 10). The rest of the lines are obtained by 
truncating the series in eq. (6.23) at linear and quartic powers of z * with a and d given in eq. 
(6.28). Similarly to the symmetry unbroken phase, figure 5, the linear approximation 2* oc l is 
better for larger d. 

Another relevant feature of the entanglement entropy eq.(6.32), that requires clarifi¬ 
cation, is that it does not obey the volume law. It is well known that for a sufficiently 
large i, the finite contribution of the entanglement entropy at finite temperature satisfies 
the volume-law, not the area, i.e., a linear-in-£ term is expected. The analytical prediction 
eq.(6.32) does not reproduce such behavior since we are neglecting terms (z*/ z 0 ) d C 1 in 
eqs.(6.23),(6.26). This is fine for small f: or, for a fixed i and T, and a sufficiently large d. 
However, for a large, but fixed d, the approximation breaks down for large £ since z* even¬ 
tually becomes sufficiently close to the horizon z 0 so that (z*/z 0 ) d ~ 1. As seen in figures 
5,6, in principle a remedy to this problem is to include more terms in the expansions given 
in eqs. (6.23), (6.26). The first correction, proportional to (z*/z 0 ) d+1 , coming from the 
n — 1 term, still leads to an analytical, but cumbersome, expression for z* (£) in the case of 
d — 3. Indeed, as is shown in figures 5,6, by including this term, the analytical expression 
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agrees with the numerical results up to larger values of l, however we do not yet observe 
the expected area law for and l —>■ oo. Indeed, for any finite number of terms the approxi¬ 
mation inevitably still breaks down at some finite i, and, already for d = 3, the subleading 
correction, oc {z* f z^) 2d ~ x , coming from the n = 1 term in eq.(6.23), leads to a fifth degree 
polynomial whose roots cannot in general be found numerically. Consequently we keep 
only the leading correction in the equations for z * and 2* so our results are fully analytical. 
As is shown in figure 7, by including this additional term only, the analytical expression 
for the finite part of the entanglement entropy agrees well with the numerical results in 
the range of i shown. However, as was mentioned previously, we do not yet observe the 
expected area law for l oo. If, on the other hand, we carry out an analogous expansion 
in the parameter 1 — (z*/z 0 ) d , see appendix B, we obtain the expected linear dependence 
of the entanglement entropy s oc i. 

We also note that the finite subleading contribution, second term of eq.(6.32), that 
does not depend on the scalar, has already been reported in Ref. [35, 42], The dependence 
on the scalar, proportional to n < 0, is consistent with previous numerical results [37]. It 
is smaller in the symmetry broken phase and its temperature dependence is not analytical 
at T c due to the different prefactors in the temperature dependence of the entanglement 
entropy in the broken and unbroken phase. We note the temperature enters both through 
zq and /i and it depends quadratically on the strip length £, for small i. These results 
illustrate the potential of 1/d expansions to obtain analytical results in problems where 
only numerical calculation were available. 
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d- 2 



Figure 7. Finite part of the entanglement entropy, s, of a rectangular strip, A, of width i in 2 
spatial dimensions (d = 2 + 1), where = *>*:?; + s and Sdi V contains the UV-cutoff. The blue 
dots are obtained by numerical integration of eqs. (6.21) and (6.22) in the Reissner-Nordstrom 
background with f(z) = /rn, zq = 0.992, /j, = 2.02, Q 2 = /x 2 Zq 7 2 , 7 = Similarly, the 

red dots, hardly indistinguishable from the blue, correspond to the integration of these equations 
in the superconducting phase where f(z) = /rn + e 2 /2 and f 2 = is given in eq. A. 10. 

We set m 2 L 2 = 0, zq = 1, q = 4 and e = 0.2 corresponding to T/T c ~ 0.995. The expansion 
parameter e = (0/(2 A — d )) 1//A is defined in the same appendix. In the superconducting phase 
/rn is given in terms of fj ,o = 2.00, eq. (A.5), and // in terms of k = —0.78, eq. (A.12). The solid 
lines are obtained from the analytical results for the superconducting and normal background from 
eqs. (6.23) and (6.26) by neglecting terms of 0(z*jzo) 2d ~ l and 0(z^/zo) 2d ~ 2 , respectively. The 
dashed lines are obtained from the linear approximation, z* = 2* = i ——/////—. As anticipated, 

2 V 7rr ( 2d—2 ) 

the analytical estimation of the entanglement entropy, calculated in the large d limit, breaks down 
as i increases. Nonetheless the qualitative behavior is very similar even for d = 2 + 1 dimensions 
in the UV-boundary. Inset: difference between the finite parts of the entanglement entropies in 
both phases as a function of the strip length. As before, the dots correspond to the numerical 
results while the dashed and continuous lines are our analytical results corresponding to the linear 
approximation and the subleading quartic correction, respectively. The lines are restricted to the 
region, in £, where the estimations of the tip of the minimal surface in each background 2*, z* < zo, 
see figures 5,6. 
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7 Conclusions 


We have studied the entanglement entropy and the conductivity in holographic supercon¬ 
ductors at zero and finite temperature in the limit of large spatial dimensionality. The 
coherence peak of the conductivity becomes narrower and the ratio between the energy 
needed to break the condensate and the critical temperature decreases as the spatial di¬ 
mensionality increases and have a well defined d —» oo limit. This is a clear indication 
that the coupling of the scalar with the bulk is weaker in the large dimensionality limit. It 
would be interesting to explore whether there is a bound for these quantities in theories 
with gravity duals. We have computed the dependence on the dimensionality d on the 
entanglement entropy at zero and close to the critical temperature and for the conductivity 
at zero temperature. Our results confirm the expectation that the entanglement entropy 
is smaller in the symmetry broken phase with a difference that decreases with the spatial 
dimensionality. These results are a strong indication that large d expansions are a helpful 
tool to obtain analytical results in holography. 
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A Appendix: Entanglement entropy at T ~ T c 

The equations of motion, expressed in the coordinate z = 1/r, are: 



(A.l) 


q 2 e x ppL A 

P 


Close to the transition, the fields can be expanded in powers of e = (3, [44, 45], where 
il)(z —> 0) ~ T( —) A , and A = A + is the larger conformal dimension. More concretely, 

^ ~ e^i + e 3 ^ 3 + ..., 0 ~ 0o + e 2 0 2 + - • •, f ^ fo + e 2 f 2 +■ ■ ■, x - ^X 2 + - ■ ■, (A.2) 
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For the purpose of the entanglement entropy calculation given in section 6.2 we calculate 
analytically the first non-trivial terms of this field-expansion in the region where ^0 j <C 
1, which, for larger d allows z to approach zo with a better level of approximation than 
for small d. However, here we compute all the terms in the perturbative expansion up to 
0{e 2 ). From the equations of 0 and / given in eq. (A.l), it is easy to see that the first 
zeroth order terms of these fields are: 


00 — d0 


1- I - 

Zo 


d-2 


fo(z) = 1 — (1 + Q 2 ) ( — ) +Q 2 ( — 


Z 0 


Zo 


2d—2 


(A.3) 


where Q 2 = do z o^ 2 an d do is the chemical potential at which the scalar field condenses. 
X = 0 and / = f 0 corresponds to the Reissner-Nordstrom black hole with planar topology. 
The equation for the first term in the expansion of 0, eq. (A.2), is well known: 

0 2 L 4 0 2 rn 2 L 2 ' 


0 = 0 " - 




/o 2 


Z 2 fo 


(A.4) 


giving the expected (z/zti) A + ... behavior close to the boundary. 0 X can be obtained by 
rewriting eq. (A.4) as a Sturm-Liouvillle eigenvalue problem, [46, 47], and using as ansatz 

, a is given by the value, a c that minimizes the following 


OLZ 


d—1 


01 = z a F 0 (z), F 0 = 1 
expression, 

fgdz z 2A ~ d+1 (l-z d ) 
M 2 (a) = - - 


771/2 

^0 


" m2L2 +A(A -d)-Ag) § 


1 — Z d 


f 0 dz z 2A ~ d+1 F 2 q 


2 ( 1 —* d ~ 2 ) 2 
l-z d 


(A.5) 


and Hq = M 2 (a c ). The equation for y 2 is: 


, 2 zi>' 2 2g 2 ^L 4 0 2 0 2 _ 


X 2 = 


d—1 


+ 


(rf-l)/o 2 


= FJz), 


(A.6) 


and thus X 2 .(z) = dz 1 F x (z r ). Similarly, from the equation for i 


d -3 


2<? 2 0O01 1 M , , _ jp f ^ 

-Jh ~ “ 2* oX2 = F * {Z) ’ 


(A.7) 


and the leading behavior of 0 2 is given by the homogeneous solution. Close to the horizon, 
02 is expected to receive corrections from the last two terms in eq. (A.7). However, 
we impose such corrections, controlled by 0 O and 0], to satisfy the boundary condition 
02(zo) = 0 independently of the homogeneous solution. Therefore, 


02 ( z ) = k 


1 - - 


02 (Z) = 02(2) +0 2 (», 

d-2" 


, 02 (z) = I duu 


"ZQ 


d- 3 


~° dv w 


(A.8) 
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Finally, the equation for / 2 is given by: 

d z 3 LV 0 < 1 

2 z 2 d — 1 z{d — 1) 


$ + 2J 2 / 0 ^f + m 2 L 2 ^l q L 


The previous equation can be integrated straightforwardly, 

f2(z) = mz)+f b 2 (z), 

PqKZqL* 
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(A.9) 


(A ' 10) 

U d 


In the largc-c/ limit, f 2 dominates over f 2 . The only parameter to be determined is k, 
which follows from the equation of wp 




d~ 1 /£ 




q z L A (f) o m 2 L 2 


V>3 = -T0O, 


r * = (| - y - ) </>; + 


m 2 L 2 
~Wo 


2 g 2 L 4 

~f3“ 

J o 


/ 2 - 


(/ 2 0O — /o0 2 ) + 


P 

Jo 


ipo- 

(A.ll) 


From the previous equation and using eq. (A.4), it follows immediately, 0 = fj° dz , 

which imposes a condition on k: 
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(A.12) 

Finally, from eq. (A.8), for some e > 0, the chemical potential is given by p ~ 

/i 0 + e 2 (/t + 0 2 (O)) + ^( e4 )’ P — P 0 + e2fi: an d the temperature T < T c : 

T f'(z 0 )e~^P 


-fo( z o) 
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d — 1 


pd- 1 4:17 (p 0 + e 2 K.) d- 1 
2 I 7 2 (-o) X 2 (-o) 


1 + e 3 
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fo(zo) 


(d- l)pofo{zo) 


(A.13) 


where, /q(z 0 ) = — d + // ( 2 ^ 2] ) . For e = 0, the previous equation gives an estimation 
for the critical temperature. This expression is more complicated that the one given in 
Ref.[16], in which the backreaction of the scalar on the geometry is neglected. Notice 
however, that, we were not after an alternate result for T c , in fact, in this section we have 
not used the largc-ri limit since we explicitly look for all the terms that modify the geometry 
close to the phase transition. In order to analytically evaluate the leading correction on the 
entanglement entropy of a strip with its complement, section 6.2, we take the leading 
correction on the blackening function, f 2 (z), given in eq. (A.10). 
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B Large t limit of the entanglement entropy at finite temperature and 
fixed d 


In the RN background, from eq. (6.22), 


1 

2 



(B.l) 


Let us split the region of integration into two: [0, z*] = [0, z a \ U (z a , 0 *], for some 0 < 
z a < 0 * and let us rename the integral in the first interval as t\ /2. In the second interval 
we change variables to z = z 0 — e and expand the integrand for e/z 0 -C 1. 
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(B.2) 


In spite of the explicit dependence of the third term on z*, in the limit z* —> zo we can take 
it as a divergent contribution, £&„, while the middle term remains finite, thus: 
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Similarly, splitting the integral of the entanglement entropy, eq. 
regions the integrations carries analogously, 


- £div (B.3) 

(6.21), into the same 


\/Zo f Zo Zif \ 1/2 sjzo f Zo - z a \ 1/2 

S ^ Sl y/(d-2)Q*-d\z2*- 2 -z* l - 2 J V(ci-2)Q 2 -Ar-^V 

(B.4) 

where si contains the UV-cutoff and the last term is also divergent in the limit z* z () . 
The middle term can be substituted using eq. (B.3) which leaves a term proportional to i. 
In the limit — » z 0 , £ — » oo however, the term in eq. (B.3) is regularized by £ div , yielding 
a finite term. 
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C Electrical conductivity at T > 0 

The boundary conditions near z —» 0 are: 
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(C.l) 

A is a cutoff which affects only the imaginary part of cr. We take A = 1. From the above 
expressions and the eq. (3.6), the conductivity is: 
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